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Frequency Offset Estimation for Communication Systems Method and Device 

for Inter Symbol Interference 



Field of the Invention 

The present invention relates to frequency offset estimations for use in 
communications systems, and particularly to a frequency offset estimation based 
on state-based differential accumulation. 

10 

BACKGROUND OF THE INVENTION 
p Many modern communication systems are carrier-modulated systems that 

% transmit information bearing signals through some specified bands. Due to the 
pf mismatch between the local oscillators of the transmitter and the receiver, a carrier 
N s 15 frequency offset is seen at the receiver. In general, a large frequency offset causes 

fil 

rapid phase rotation, which could be far beyond the tracking capability of the 
5 receiver functions, such as the carrier phase tracking loop and the channel 

jj .j \ 

% estimator. It is thus desired to have an initial frequency offset estimation to reduce 

O the amount of frequency offset seen by the remaining part of the receiver. Some 
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20 current frequency offset estimation algorithms are based on the ideal channel 
condition of additive white Gaussian noise (AWGN) and in some applications, a 
communications channel experiences severe inter-symbol interference (ISI). One 
example is the multi-path fading channel for wireless communications where ISI 
leads to extra noise terms and degrades the performance of prior art frequency 

25 offset estimations. 

Several prior art frequency offset estimation schemes used in the ISI channels 
are briefly described herein. A "Fourier Coefficients after Non-linear Operation" 
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algorithm performs the following steps: 

■ Do some non-linear operation, such as differential detection, 

■ Calculate the Fourier coefficients for some fundamental frequencies at the 
output of the non-linear operation, 

■ Perform another operation on these coefficients. 

■ Estimate the frequency offset from its phase. 

This algorithm does not require knowledge of received data symbols and timing, 
but it makes some approximations based on the assumption of the channel 
response. Another prior art scheme, "Jointly Channel and Frequency Offset 
Estimation", utilizes a maximal-likelihood (ML) type of jointly channel and 
frequency offset estimator in which the algorithm has to search the maximal point 
of the non-linear ML equation after certain matrix operations are performed. This 
algorithm requires the knowledge of received data symbol and timing. The "Phase 
Rotation of Channel Estimation" method estimates frequency offset by neglecting 
the effect of frequency offset within a short time, then the frequency offset is 
estimated by the phase rotation of the channel coefficients. This method splits the 
channel and frequency offset estimations and requires knowledge of received data 
symbols and timing. The "Maximum Stated-based Accumulation (MSA)" method 
utilizes a modified version of the non-linear ML equation which is obtained by 
collecting the terms corresponding to the same state to form a whole while finding 
the maximal point of the modified non-linear equation. Knowledge about the 
received data symbols and timing is required. 

SUMMARY OF THE INVENTION 
The object of the present invention is to provide a frequency offset estimation 
system and method utilizing inter-symbol interference in a channel 
communication system. An advantage of the present invention is that it does not 
require intensive computation of Fourier coefficients or maximization of any non- 
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linear equations as in the prior art. Further, according to the present invention, it is 
not necessary to frequency offset-estimate the channel jointly or concurrently 
while the differential operation over the received samples of the same state is 
being performed. Thus the present invention is simpler and easier to implement 
5 than other prior art techniques. 

The present invention provides a frequency offset estimation system and 
method utilizing a slow channel variation with respect to a predetermined 
observation time, the channel being treated as in a static condition within an 
observation window. In a receiving unit of the channel, there are a plurality of 
10 data symbols being pre-stored or detected and a receiving sample sequence being 
received. A predetermined number of data symbols of the plurality of data 
% symbols are defined as a state, and the plurality of data symbols are divided into a 
W first, a second, and up to an / state, each having a length equal to a 
fl| predetermined number whereby that predetermined number is smaller then the 

U 15 number of said data symbols. The receiving symbol sequence has a phase shift 
and an angular frequency offset and is formed by a plurality of receiving sample 
2 sections. Each of said plurality of receiving sample sections contains a plurality of 

W samples which are received continuously. The method of the present invention 
0 comprises the steps of first executing a complex multiplication of two selected 

20 receiving sample sections of the plurality of receiving sample sections of the 
receiving sample sequence so at to obtain a value wherein the two selective 
receiving sample sections are mapped to the same state of the data symbol and 
each receiving sample section is designated by a symbol index whereby the 
symbol indices form a series of continuous integers corresponding to the 
25 respective continuous receiving samples of the received symbol sequence; then 
dividing the multiplied value by a difference of the respective symbol indices of 
the two selective receiving sample sections thereby acquiring an estimated 
frequency offset value; and finally weight-averaging a plurality of estimated 
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frequency offset values for different receiving sample section pairs. 

The various objects and advantages of the present invention will be more 
readily understood from the following detailed description when read in 
conjunction with the appended drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 

Fig. 1 shows the hardware implementations and block diagrams for the SDA 
algorithm without state reduction according to the present invention. 

Fig. 2 shows the hardware implementations and block diagrams for the SDA 
algorithm with state reduction according to the present invention. 

DETAILED DESCRIPTION OF THE INVENTION 

The preferred algorithmic embodiment of the present invention is termed the 
"State-based Differential Accumulation (SDA)" and utilizes knowledge of 
received data symbols and timing. In the following description of the present 
invention, the concept of the state-based operation is first described by 
mathematical derivation. Then the detailed operations are presented utilizing both 
equations and conceptual descriptions. 

Assuming that the channel variation in a communications system is slow with 
respect to an observation time, and a channel is thus treated as in static condition 
within an observation window, then a continuous-time received signal can be 
written as, 

r{t) = A e jff e J ^ *(' " kT "> + w & 

k=~O0 

where 

A is a scalar and is the gain of the channel, 
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6 is the phase shift of the received signal from the transmitted signal, 
co d is the angular frequency offset of the channel, 
d k is the k-th transmitted data symbol. 

g(t) is the equivalent composite impulse response combing the effects of 
5 transmission pulse shaping filters, channel impulse responses, and receiving 
filters. 

T is the period of data symbol transmission. 
w{i) is the additive noise process. 

Taking samples at a rate of ^~ , a received sample sequence is obtained 

T 

\* 10 r[n) = r(t = n-—) = Ae J *e J " dn "fd k g(n' — -kT) + w(t = n-—) 

Q ' M {T K M M 



|j =Ae je e " £^(_.(„-H/)) + h</ = ».— ) 



./•— * 



A" =-oo 

where 

= <z; d • T is the angular shift over one symbol period, and 



fU 15 g[n] = g(t = n~) 

M 

T 

^n] = w(t = n-—) 
M 

where g[n] and w[n] are the discrete-time sampled version of their continuous- 
time counterparts. 

In the concept of state-based differential operation, for practicality, the 
20 impulse response is truncated to some finite length, for example L times the 
symbol periods. Without loss of generality, assuming that 
gM = 0 Vne {0,1,2,..., (LAf-1)} 
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then consider « = / • M + i , where leZJe {0,1,2,..., (M - 1)} , 

i 

*=/-£ 



r[« = M + i] = A e Jt e Ja ' (f * M) £ rf t - £)M + i] + 



/n </ (/+— ) 
= e " -jr(S /f i) + w[n] 



5 where 



jr(5 l ,o=^2;^g[(/-*)w+i] 



k=l-L 



s i = {^/-i > ^/-i+i , • ■ • , ^ } is the state during the / -th symbol. 



M< Then the state is formed by the current and previous L data symbols. Only these 

b 

rt 10 symbols contribute to the channel ISI related term, X(S l9 i). The receiving 



samples of the same state and the same offset / can be used to remove the channel 
effect and estimate the frequency offset. This can be shown as follows. 
Neglecting the noise term, for any S h = S h and l x * l 2 , gives 

X(S h J) = X(S l2 ,i). 

15 Performing the differential operation and defining a new variable for the result, 
gives 

utfM = r *[l x M +i]r[l 2 M +i] = e* M) . | ,/) | 2 

The angle of the differential operation result contains the frequency offset 
information. Thus, the frequency offset can be obtained by 

20 n^-LzutVM 

A/ 

where A, = / 2 - l x is the index difference. 

To reduce the variance of the estimator, multiple occurrences at the output 
of the differential operator can be accumulated before taking the phase. For 
different values of the index difference A,, the accumulation is preferably done 
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separately since they are not coherent in phase. The accumulation variable can be 
calculated as 

M-l 

v ^i)= Z Z "MM 

t 2 ~li=Ai r=0 

For any two indices spaced with A, and having the same state, a differential 
operation is executed with their corresponding samples and the results are 
accumulated according to the accumulation variable of the index difference A, . 
Frequency offset can be estimated from each of the accumulation variables by 

A / 

Due to the ambiguity from phase wrap-around, the resolvable range of the 
estimator is 

l^l<f • 

This leads to a consideration to avoid using large values of A, . 

A weighted average can be performed over the frequency estimation from 
different values of A n 

£C(A,) ^ 

where 

C(A ; ), A, = 1,2,..., A /max are the weighting factors. 

A maximal value of A, is set under the consideration of frequency estimation 
range, as described above. The selection of weighting factor preferably minimizes 
the variance of estimation, and it may also take implementation complexity into 
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account. In a simplest way, equal weights can be applied. Alternatively, the 
weighting factor can be proportional to the relative frequency of occurrences of 
A, . If the statistics of the data symbol sequence is known and fixed, the weighting 
factors can be calculated in advance. Otherwise, a set of counters is used to 
dynamically count the number of the occurrences of each A, . Then the counter 

values are used to generate the weighting factors. 

The method of the present invention can be summarized as the following rules; 

■ Execute a complex multiplication of two segments in a receiving sample 
sequence, wherein the two segments are mapped to the same state of a 
series of data symbols and each segment is indicated by a symbol index; 
and the symbol indexes are a series of continuous integers corresponding 
to respective continuous segments of the received symbol sequence; 

■ Divide the multiplied value by a difference value of respective symbol 
indexes of the two segments so as to acquire an estimated frequency offset 
value; and 

Obtain an average estimated frequency offset value by weighted-averaging a 
plurality of estimate frequency offset values for different segments. 

The number of the states is M d L+ \ where M d is the size of the set of the data 

symbol. The present invention preferably reduces the number of states. Consider 
a rotated received samples by the complex conjugate of the data symbol, and 
neglecting the noise term; 

r R [n = IM + i] = d* ♦ r[n = IM + i] 

= A e J& e aAl+ ^ £ d*d k g[(l - k)M + i] 

k=l-L 

jQ d (/+— ) 

= e » -X(S RJ ,i) 

where 
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S Rl ={d l *d i _ L9 d g *d l _ L¥l9 — 9 d l d l } = d l •{d l _ L9 d l _ M9 —,d i } = di -S t 
It is seen that phase rotation to the received samples causes the same phase 
rotation to the states. Also note that the state after rotation becomes 
^ R) i={d*d l _ L9 d*d l _ L+l9 -' 9 d l </,} 

= {di*d l _ L9 d*d i _ M9 — 9 \d l | 2 }. 

In the case of PSK modulation, the last element in the state is a constant. The 
number of the states becomes M d L , and the reduction factor is M d . In a more 
general way, the rotation phase can also be a mapping function from the data 
symbol other than its complex conjugate. 

For the case of spread spectrum, if the spreading code or spreading waveform 
is periodic with the same period as the symbol time, the spreading waveform can 
be viewed as part of the pulse shaping and can be absorbed into the equivalent 
composite impulse response. Since the method of the present invention does not 
require the knowledge of pulse shaping or channel impulse response, it is also 
applicable to spread spectrum communication systems. 

A block diagram without state reduction is shown in Figure 1. The operations 
in the block diagrams are explained as follows. The detailed and mathematical 
descriptions of each operation are described above. Fig. 1 illustrates a frequency 
offset estimation circuit according to the present invention. The circuit comprises 
a state selector 102, a receiving end 103, a state bank 107, a differential operator 
111 and an accumulating and averaging unit 114. The state selector 102 
determines a current state as a selective state (1011) based on the current and some 
previous data symbols (101) which is pre-stored in the receiving end 103 or is 
roughly detected by the receiving end 103, the state selector outputs the selected 
state (1011). The receiving end 103 serves for receiving a receiving sample 
sequence (104) and having a counter 105 for generating a symbol index (106) 
responsive to the input sample sequence (104). The receiving end outputs the 
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received sample sequence (104) and the symbol index (106). The state bank 107 
has inputs of the selected states (1011) from the state selector 102, received sample 
sequence (104) from the input end 103 and a symbol index (106) from the counter 
105. The state bank 107 divides the received sample sequence (104) into a 

5 plurality of data sample sections responsive to the selected state (1011). Each state 
unit stores one of the plurality of sample sections and a symbol index of last-time 
visit. Only contents of a selected state are modified responsive to the selected 
state (1011) and the symbol index (106). Then a selected sample section and a 
selected symbol index of the last-time visit are shifted to an output end of the state 

10 bank 107. The contents of the state are then updated by new values of the sample 
sections and the symbol index. The state bank 107 outputs the samples sections 
and the selected symbol index of the last-time visit. The differential calculator 111 
is installed between the state bank 107 and the receiving end 103 and receives 
inputs of the received sample sequence (104) and the symbol index (106) from the 

15 receiving end 103 and the selected sample section and the symbol index from the 
state bank 107. The differential calculator 111 determines a complex conjugate to 
the selected samples of the selected last-time visit through a conjugate calculator 
113. Then it executes a complex multiplication to the received samples and the 
selected sample section of last-time visit to acquire a multiplied value through a 

20 complex multiplier 112. The difference between the symbol index from the 
receiving symbol sequence and the symbol index from the state bank is also 
calculated to obtain an index difference. The differential calculator outputs the 
multiplied value and the index difference. The accumulating and averaging unit 
114 has inputs of the multiplied value and the index difference from the 

25 differential calculator 111. The accumulating and averaging unit 114 has a 
plurality of branches each corresponding to one index difference. The inputted 
multiplied value is dispatched to a corresponding branch 115 responsive to the 
index difference. The phase angle of the multiplied value is calculated by, for 
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example, an inverse-tangent operation. The phase angles are scaled differently. 
The phase angle of each branch is divided by the index difference value of the 
branch. Finally, a weighted average operation is performed over the branches by a 
weighted average unit 119, thereby acquiring the frequency offset in the receiving 
unit. The state bank 107 further comprises a distributor 108 (in the preferred 
embodiment, a multiplexer), a plurality of state units 1071 and a simple selector 
110 (in the preferred embodiment, a demultiplexer). The distributor 108 has 
inputs to the received sample sequence 104 and the symbol index 106 from the 
receiving end 103. It divides the sample sequence 104 into a plurality of sample 
section responsive to said symbol index 106 for distributing the sample sections 
and symbol indexes 106 to a corresponding state unit 1071 responsive to the 
selected state 1011 from the state selector 102. The plurality of state units receives 
and stores corresponding inputs of the data sample sections and the symbol index. 
Each state unit 1071 has a symbol memory 1072 for recording the sample section 
and an index memory 1073 for storing the symbol index 106 from the receiving 
end in response to the selected state. The symbol selector 110 serves for selecting 
the data symbol and the symbol index of the last-time visit responsive to the 
selected state from the state selection. Each of the branch further comprises an 
accumulator 116, a phase calculator 117 and a scale calculator 118. The 
accumulator 116 has an input of the multiplied value for accumulating the 
multiplied values and the index difference. The phase calculator serves for 
determining a phase of the multiplied value. The scale calculator 118 serves for 
determining a scale of the multiplied value. 

A block diagram with state reduction is shown in Figure 2. A state reduction 
phase means 120 has an input of a selected data symbol from the plurality of data 
symbols 101 for calculating a phase of the selected data symbol, and outputs the 
phase of the selected data symbol. A first phase rotation means 121 receives the 
phase of the selected data symbol for rotating phases of other data symbols with 
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the phase of the selected data symbol. A second phase rotation means 122 
receives the phase of the selected data symbol for rotating phases of the receiving 
symbol sequence with the phase of the selected data symbol. As described above, 
the complex conjugate of the data symbol can be selected as the state-reduction 
phase. Other predetermined mapping from the data symbol to the phase is also 
possible. In general, the phase rotation block can be performed by the complex 
multiplication operation. When the phase rotation is only 0 or 180 degrees, only 
the polarity inverse is needed for phase rotation. When the phase rotation is 90 or 
270 degrees, the real and imaginary parts are exchanged to perform the polarity 
inverse. 

Although the present invention has been described with reference to the 
preferred embodiments, it will be understood that the invention is not limited to 
the details described thereof. Various substitutions and modifications have been 
suggested in the foregoing description, and others will occur to those of ordinary 
skill in the art. Therefore, all such substitutions and modifications are intended to 
be embraced within the scope of the invention as defined in the appended claims. 
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